DISCOVERY OF SPECTRAL TRANSITIONS FROM TWO 
ULTRA-LUMINOUS COMPACT X-RAY SOURCES IN IC 342 



A. Kubota/ T. Mizuno,^ K. Makishima/ Y. Fukazawa,^ 
J. KOTOKu/ T. Ohnishi/ and M. Tashiro^ 
l:Department of Physics, University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo, Japan 113-0033 

e-mail : aya@amalthea.phys.s.u-tokyo. ac.jp 
2:Department of Physical Science, Hiroshima University, 1-3-1 Kagamiyama, 
Higashi-Hiroshima, Japan 739-8526 



Received 



accepted 



- 2 - 



ABSTRACT 



Two ASCA observations were made of two ultra- luminous compact X-ray sources 
(ULXs), Source 1 and Source 2, in the spiral galaxy IC 342. In the 1993 observation, 
Source 2 showed a 0.5-10 keV luminosity of 6 x 10'^^ ergs s^"*^ (assuming a distance 
of 4.0 Mpc), and a hard power-law spectrum of photon index ~ 1.4. As already 
reported, Source 1 was ~ 3 times brighter on that occasion, and exhibited a soft 
spectrum represented by a multi-color disk model of inner-disk temperature ~ 1.8 
keV. The second observation made in February 2000 revealed that Source 1 had made 
a transition into a hard spectral state, while Source 2 into a soft spectral state. The 
ULXs are therefore inferred to exhibit two distinct spectral states, and sometimes 
make transitions between them. These results significantly reinforce the scenario which 
describes ULXs as mass-accreting black holes. 



Subject headings: black hole physics — galaxies: spiral — X-rays: galaxies 
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1. Introduction 

Arm regions of nearby spiral galaxies have long been known to harbour luminous point-like 
X-ray sources (Fabbiano 1989, 1998; Read et al. 1997; Roberts & Warwick 2000). Except for 
several identified with young supernova remnants, they remained unidentified in other frequencies. 
Being often time variable, these sources are thought to be accreting single objects; we hence call 
them ultra- luminous compact X-ray sources (ULXs; Makishima et al. 2000, hereafter Paper I). 
With the X-ray luminosity reaching 10^^~^° ergs s~^, the ULXs are suspected to contain stellar 
black holes (BHs) of mass up to ~ 100 Mq or more (e.g., Colbert et al. 1995; Fabbiano 1998; 
Immler et al. 1999; Roberts &; Warwick 2000), so as not to violate the Eddington limit. However, 
this much exceeds the mass estimates (5 — 15 M©; Tanaka & Lewin 1995) for accreting stellar BHs 
in the Galaxy and the Large Magellanic Cloud. There is no consensus on the presence of ~ 100 Mq 
BHs, even though such medium-mass BHs have been proposed under slightly different context 
(Colbert & Mushotzky 1999; Matsumoto & Tsuru 1999; Ptak & GrifHths 1999). Furthermore, 
we do not find ULXs in the Milky Way or M31. Therefore, the BH scenario of ULXs remained 
inconclusive, and their nature has remained a big mystery. 

A great leap has been achieved through spectroscopy of a dozen ULXs with ASCA (Takano 
et al. 1994; Okada et al. 1998; Colbert & Mushotzky 1999; Mizuno et al. 1999; Mizuno 2000; 
Kotoku et al. 2000; Paper I), as the spectra of majority of them have been reproduced by 
so-called multi-color disk model (MCD model; Mitsuda et al. 1984) that describes emission 
from an optically-thick standard accretion disk around a BH (Shakura & Sunyaev 1973). These 
results significantly reinforce the BH scenario of ULXs, since the MCD model can explain the soft 
X-ray spectra of Galactic/Magellanic black-hole binaries (BHBs) in "soft (or high) state" (e.g., 
Makishima et al. 1986; Ebisawa et al. 1993; Dotani et al. 1997; Kubota et al. 1998; Feroci et 1. 
1999; Paper I). 

The scenario is nevertheless still clouded by a self-inconsistency, that the measured inner-disk 
temperature of ULXs, Tin = 1-0 — 1.8 keV, is too high for the implied high BH mass (Colbert & 
Mushotzky 1999; Paper I). Although the scenario may be salvaged by invoking rapid BH rotation 
(Mizuno et al. 1999; Paper I; Mizuno 2000; also Zhang et al. 1997b) and/or "slim disk" concept 
(Abramowicz et al. 1988; Watarai et al. 2000), it is important to assemble further observational 
clues. 

Here, we report further evidence supporting the BH interpretation of ULXs, i.e., transitions 
of some ULXs between two distinct spectral states. These results are based on two observations 
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with ASCA (Tanaka et al. 1994) of two ULXs, both in the spiral galaxy IC 342 at a distance of 
4.0 Mpc (Tully 1988). 

EDITOR: PLACE FIGURE |l| HERE. 



2. Observation 

The first observation of IC 342 was performed on 1993 September 19, as reported by Okada 
et al. (1998). The obtained X-ray image is reproduced in Figure |l]a. In addition to unresolved 
emission from the nuclear region, we see two ULXs, the brighter Source 1 and the fainter Source 2, 
both known from the Einstein era (Fabbiano & Trinchieri 1987). On this occasion, Source 1 
exhibited clear short-term variability (Okada et al. 1998). 

We observed IC 342 with ASCA again on 2000 February 24-March 1. The CIS was operated 
in the normal PH mode, while the SIS in the Faint mode. Although the acquired X-ray image. 
Figure ||b, again reveals the two sources, their relative intensities have reversed from those in 
1993. As quantified in Table 1, the Source 1 flux decreased to ~ 30% of its 1993 value, while that 
of Source 2 increased by a factor of 1.8 meantime. 

EDITOR: PLACE TABLE | HERE. 



3. Data analysis and results 
3.1. Data in 1993 

Spectral results on Source 1 from the 1993 observation are already described in Okada et al. 
(1998), Paper I, and Mizuno (2000). Referring to Paper I, we simply quote here that the 1993 
spectra of Source 1 can be reproduced successfully by an absorbed MCD model of Tin ~ 1.8 keV. 
In Table 1, we reproduce essence of these results. 

We here analyze the GIS/SIS spectra of Source 2 after Mizuno (2000). We selected good data 
in the same manner as in Okada et al. (1998), and accumulated photons around the Source 2 
centroid, separately for the four detectors. We subtracted background using blank-sky data. We 
then added the SISO and SISl spectra into a single SIS spectrum, and those from GIS2 and GIS3 
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into a single GIS spectrum. The net exposure became 38 ks for the GIS and 36 ks for the SIS. 
The derived Source 2 spectra, shown in Figure ||a, are much harder than those of Source 1, and a 
joint GIS/SIS fit with an absorbed power-law model of photon index T ~ 1.4 has been successful 
(Table 1). Although the absorbed MCD fit is also acceptable, it requires an unrealistically high 
disk temperature of Tin ~ 3.0 keV (Table 1). An absorbed Bremsstrahlung model of temperature 
> 23 keV is also acceptable {x^/^ = 101.1/87). Considering typical spectra of BHBs, we consider 
the power-law fit to be most appropriate for Source 2. 

EDITOR: PLACE FIGURE | HERE. 
3.2. Data in 2000 

We use only the GIS data for our spectral study of the 2000 observation, because the 
long-term SIS degradation has made its low energy response rather uncertain. We selected good 
GIS data using criteria of geomagnetic cutoff rigidity > 6 GeV, and the target elevation above the 
earth rim > 5°. This has yielded 276 ks of good exposure. 

Through the long observation. Source 1 and Source 2 both showed a relatively steady 0.7-10 
keV intensity of 0.027 c s"^ and 0.045 c s"^ per GIS detector respectively, with mild variation. 
We here utilize the whole length of GIS data for spectral evaluation, by accumulating events over 
two circular regions of radius 2. '5 around Source 1 and 3.'0 around Source 2. By subtracting 
background utilizing blank-sky data, and adding data from the two GIS detectors, we have 
obtained the Source 1 and Source 2 spectra, shown in Figures ^3 and |2|c, respectively. 

The Source 1 spectrum has become much harder and less convex than in 1993, and fitted 
approximately {x^ = 119.8/88) with an absorbed power-law of F 1.8. However at 7-10 keV, 
we observe negative residuals suggestive of an Fe-K edge. By introducing an edge absorption 
at ~ 8.4 keV, the fit became acceptable (Table 1); an F-test indicates that the edge feature is 
significant at 99.5% confidence. Although the MCD fit is not too bad (Table 1), the absorption it 
requires falls below the Galactic value of ~ 3 x 10^^ cm~^. We therefore consider the power-law 
fit with an Fe-K edge appropriate. 

In contrast, the Source 2 spectrum has become significantly more convex than in 1993. It 
can now be expressed adequately with an MCD model of Tin = 1-6 keV, whereas the power-law 
fit is unacceptable (Table 1). Thus, the Source 2 spectrum in 2000 is close in shape to the 
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Source 1 spectrum in 1993, except for a higher photoelectric absorption. An inclusion of a separate 
power-law component, with photon index fixed at 2.5, yields only an upper limit of 15% in terms 
of its 0.5-10 keV flux, and the MCD parameters are not affected beyond the fitting errors (see also 
§ 3.3 of Paper I). This indicates dominance of the optically-thick disk emission in the ASCA band. 

To examine possible changes in the absorption, we fitted the Source 1 spectra in 1993 and that 
in 2000 with the MCD model and the power-law model respectively, under a common absorption. 
We obtained an acceptable {Xu — 1-13) joint fit, with Nn = 5.3 ib 0.5 cm~^ and little changes in 
the other parameters. Similarly, the 1993 power-law fit and the 2000 MCD fit to the Source 2 
spectra accepted {Xu — 0.94) a common absorption of A'^h = 18.1 it 0.8 cm^^. Therefore, the 
transition of neither source was accompanied by measurable changes in the absorption. 

4. Discussion 

In the first ASCA observation of IC 342, Source 1 exhibited a high luminosity exceeding 
-|^q40 gj,gg g-i (^'Pable 1), prominent short-term variation, and a convex spectrum expressed with 
an MCD model of a high disk temperature. Tin ~ 1.8 keV. In Paper 1, these features have been 
identified as prototypical ULX properties. 

In the same observation. Source 2 showed a distinct hard spectrum expressed with a 
power-law of P 1.4. Among the ULXs so far studied with ASCA, another object, i.e.. Source A 
in NGC 1313 observed in 1993, exhibited a similarly hard spectrum with P ~ 1.7 (Petre et al. 
1994; Colbert & Mushotzky 1999; Mizuno 2000); the remaining objects showed the MCD-type 
spectra (Paper I). We hence presume that there are two subtypes of ULXs, a majority exhibiting 
the soft MCD-type spectra, and a minority showing the hard power-law spectra. However, from 
these results alone, we cannot tell whether or not these subtypes represent intrinsically distinct 
two source populations. 

The second ASCA observation made in 2000 revealed that Source 1 had made a soft-to-hard 
state transition, whereas Source 2 a hard-to-soft one, although their reciprocal behavior must be 
a chance coincidence. To visualize this, we present in Figure ^ their deconvolved spectra. These 
results immediately imply that; (1) the two ULXs must be accreting compact objects, to make 
such drastic spectral changes in 7 years; (2) the two ULX subtypes represent two well-defined 
states (soft state and hard state) of the same source population; (3) some (if not all) ULXs 
make transitions between the two spectral states; and (4) at least in the ASCA band, the source 
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is less luminous while in the hard state. These properties give a conclusive support to the BH 
interpretation of ULXs, because they are just what characterize the known BHBs (e.g., Maejima 
et al. 1984; Tanaka & Lewin 1995; Tanaka & Shibazaki 1996; Zhang et al. 1997a; Wilms et 
al. 2000). A similar (though less convincing) example may be the aforementioned Source A in 
NGC 1313, because its spectrum turned very soft in 1995 (F ~ 2.8; Colbert & Mushotzky 1999; 
Mizuno 2000). 



EDITOR: PLACE FIGURE | HERE. 



Thus, the long-lasting puzzle of ULXs are being unveiled, but a series of new questions are 
arising in turn. Which is the right solution to the "too high Tin" problem? Are the inferred 
high BH masses of ULXs consistent with the current understanding of stellar BH formation? 
Why are ULXs absent in Our Galaxy and M31? Do they form a population distinct from the 
Galactic/Magellanic BHBs? In short, the ULX formation scenario is yet to be clarified. 

The authors thank the ASCA team members. The present work is supported in part by 
the Grant-in-Aid for Center-of-Excellence, No. 07CE2002, from Ministry of Education, Science, 
Sports and Culture of Japan. 
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Table 1: The spectral parameters of IC 342 Source 1 and Source 2, with 90% confidence limits. 



Epoch 


j,a) 
Jx 




Fit with the MCD model 

Tin (keV) x^jv 


Fit with the power-law model 
Photon index N'^ ")d 


Source 1 


















1993 


10.2 


1.9 


1.77 ±0.05 


4.7 ±0.3 


137.4/135 


1.90 ± 0.05 


9.3 ±0.6 


266.5/135 


2000 


3.1 


0.6 


2.06 ± 0.08 


1.9 ±0.4 


123.9/89 


1.73 ± 0.06^) 


6.4 ± 0.7^) 


101.1/86^) 


Source 2 
1993 
2000 


4.1 
7.2 


0.8 
1.4 


3.03 ± 0.30 
1.62 ± 0.04 


9.9 ±0.9 
18.2 ±0.8 


96.4/87 
90.5/89 


1.39 ±0.10 
2.48-'') 


14.3 ± 1.6 
31.8-^) 


102.8/87 
232.0/89 



°) The 0.5-10 keV source flux at the top of atmosphere, in unit of 10 ergs s ^ cm ^. 

''^ The absorption-uncorrected 0.5-10 keV luminosity in 10^° ergs s~^, assuming 4.0 Mpc 

distance and an isotropic emission. 
'^^ Column density for absorption assuming solar abundances, in unit of lO^-*^ cm~^. 

Results taken from Paper I. 
^) An ionized Fe-K edge (at 8.4 ± 0.3 keV) is apphed, with an optical depth of 0.9 ± 0.5. 

Errors are not shown since the fit is highly unacceptable. 



This manuscript was prepared with the AAS macros v4.0. 




Fig. 1. — The 0.7-10 keV X-ray images of IC 342 taken with the ASCA GIS, shown after smoothing 
with a Gaussian kernel of cr = 0'.5 and superposed on the optical image. The background has not 
been subtracted, and the X-ray contour levels are linear. The source to the lower right of the 
optical nucleus is Source 1, while that to the upper right is Source 2. (a) The data taken on 1993 
September, (b) Those of 2000 February. 




Energy (keV) 

Fig. 2. — The ASCA spectra of IC 342 Source 1 and Source 2, together with prediction of the 
best-fit model (Table 1). The model parameters are given in Table 1. Fit residuals are shown for 
both the power-law and MCD models, (a) The SIS (open circles) and GIS (filled circles) spectra 
of Source 2 in 1993, jointly fitted with an absorbed power-law model, (b) The GIS spectrum of 



- 11 - 



Source 1 in 2000, fitted with an absorbed power-law model incorporating the ionized Fe-K edge, 
(c) The GIS spectrum of Source 2 obtained in 2000, fitted with an absorbed MOD model. 




Energy (keV) Energy (keV) 



Fig. 3. — The ASCA GIS spectra of IC 342 Source 1 (panel a) and Source 2 (panel b), presented 
after removing the instrumental responses and the photoelectric absorption. Open and filled circles 
indicate the data taken in 1993 and 2000, respectively. The Source 1 spectrum in 1993 and that of 
Source 2 in 2000 are deconvolved using the MCD-model fits, whereas the others using the power-law 
fits. The solid lines indicate these best-fit models in their incident forms. 



